Can hadronic rescattering explain the large elliptic flow and small HBT radii 
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Results from the data obtained in the first physics run of the Relativistic Heavy Ion Collider 
(RHIC) have shown suprisingly large elliptic flow and suprisingly small HBT radii. Attempts to 
explain both results in a consistant picture have so far been unsuccessful. The present work shows 
that a thermal model + hadronic rescattering calculation can explain both elliptic flow and HBT 
results from RHIC. The calculation requires a very early hadronization time of about 1 fm/c after 
the initial collision of the nuclei. 
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PACS numbers: 25.75.Gz, 25.75.Ld 

Results of the Year-1 running of the Relativistic Heavy 
Ion Collider (RHIC) for Au+Au collisions at ^/s = 
130 GeV have shown suprisingly large pion elliptic flow 
and suprisingly small radii from two-pion Hanbury- 
Brown-Twiss interferometry (HBT) ||. Attempts to 
explain both results in a consistant picture have so far 
been unsuccessful. Hydrodynamical models agree with 
the large elliptic flow seen in the RHIC data [Q but sig- 
nificantly disagree with the experimental HBT radii 
On the other hand, relativistic quantum molecular dy- 
namics calculations which include hadronic rescattering, 
for example RQMD v2.4 j^, significantly underpredict 
the elliptic fiow seen in the RHIC data |7] but predict 
pion HBT radii comparable to the data A calcula- 
tion has recently been made to extract HBT radii with 
a hydrodynamical model coupled with a hadronic rescat- 
tering afterburner with the result that the HBT radii are 
significanly larger than measurements This lack of 
a single model to explain both results has been our first 
big mystery from RHIC. It has been suggested that we 
should call into question our current understanding of 
what information pion HBT measurements give us pO[ . 

In an effort to address this mystery, the present work 
explores a somewhat different picture of the nuclear 
collision than those presented above. In this picture, 
hadronization into a thermal equilibrium state occurs 
soon (about 1 fm/c) after the initial collision of the nuclei 
followed by hadronic rescattering until freezeout. The 
goal will thus be to test whether hadronic rescattering 
alone can generate enough general fiow in the system 
to explain both the elliptic fiow and HBT results from 
RHIC. Note that this approach, while similar to the rel- 
ativistic quantum molecular dynamics calculations men- 
tioned above, differs from them in the choice of the model 
of the initial state of the system before hadronic rescat- 
tering commences (e.g. the other models use a color 
string picture for the initial state) The price for us- 
ing the present method is that hadronic-like objects must 
exist during the high-energy-density (p > 1 GeV/ fm^) 
phase of the collision. The discussion of this important 
point is deferred until later. 



A brief description of the rescattering model calcula- 
tional method is given below. The method used is simi- 
lar to that used in previous calculations for lower CERN 
Super Proton Synchrotron (SPS) energies jllj. Rescat- 
tering is simulated with a semi-classical Monte Carlo 
calculation which assumes strong binary collisions be- 
tween hadrons. The Monte Carlo calculation is carried 
out in three stages: 1) initialization and hadronization, 
2) rescattering and freeze out, and 3) calculation of ex- 
perimental observables. Relativistic kinematics is used 
throughout. All calculations are made to simulate RHIC- 
energy Au+Au collisions in order to compare with the 
results of the Year-1 RHIC data. 

The hadronization model employs simple parameteri- 
zations to describe the initial momenta and space-time 
of the hadrons similar to that used by Herrmann and 
Bertsch ||l2|. The initial momenta are assumed to follow 
a thermal transverse (perpendicular to the beam direc- 
tion) momentum distribution for all particles, 

{l/mT)dN/dmT = CmT/[exp (mr/T) ± 1] (1) 



where rriT = \/pT^ + wo^ is the transverse mass, pt is 
the transverse momentum, toq is the particle rest mass, C 
is a normalization constant, and T is the initial temper- 
ature of the system, and a gaussian rapidity distribution 
for mesons, 



dN/dy ^Dexp [-{y - yof/{2ay^)] 



(2) 



where y = 0.5 In [{E + pz)/{E — p^)] is the rapidity, E is 
the particle energy, Pz is the longitudinal (along the beam 
direction) momentum, D is a normalization constant, yo 
is the central rapidity value (mid-rapidity), and ay is the 
rapidity width. Two rapidity distributions for baryons 
have been tried: 1) fiat and then falling off near beam 
rapidity and 2) peaked at central rapidity and falling 
off until beam rapidity. Both baryon distributions give 
about the same results. The initial space-time of the 
hadrons for 6 = fm (i.e. zero impact parameter or 
central collisions) is parameterized as having cylindrical 
symmetry with respect to the beam axis. The transverse 



particle density dependence is assumed to be that of a 
projected uniform sphere of radius equal to the projec- 
tile radius, R {R ^ tqA^/"^, where tq = 1.12 fm and A 
is the atomic mass number of the projectile). For 6 > 
(non-central collisions) the transverse particle density is 
that of overlapping projected spheres. The longitudinal 
particle hadronization position (zhad) and time {tuad) are 
determined by the relativistic equations [Esl, 
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Zhad = Thad sinhy; t^ad = Thad cosh 2/ 



(3) 



where y is the particle rapidity and Thad is the hadroniza- 
tion proper time. Thus, apart from particle multiplici- 
ties, the hadronization model has three free parameters 
to extract from experiment: ay, T and Thad- The hadrons 
included in the calculation are pions, kaons, nucleons and 
lambdas (tt, K, N, and A), and the p, w, ?7, 7]' , 0, A, and 
K* resonances. For simplicity, the calculation is isospin 
averaged (e.g. no distinction is made among a 7r+, tt*^, 
and TT"). Resonances are present at hadronization and 
also can be produced as a result of rescattering. Initial 
resonance multiplicity fractions are taken from Herrmann 
and Bertsch |2|, who extracted results from the HELIOS 
experiment |1J|. The initial resonance fractions used in 
the present calculations are: t\I-k = 0.05, p/tt = 0.1, 
pjuj = 3, 4>l{p + tj) = 0.12, r( It] = K*/uj = 1 and, for 
simplicity, /S./N = 0. 

The second stage in the calculation is rescattering 
which finishes with the freeze out and decay of all par- 
ticles. Starting from the initial stage {t = Q fm/c), the 
positions of all particles are allowed to evolve in time in 
small time steps {dt = 0.1 fm/c) according to their ini- 
tial momenta. At each time step each particle is checked 
to see a) if it decays, and b) if it is sufficiently close to 
another particle to scatter with it. Isospin- averaged s- 
wave and p-wave cross sections for meson scattering are 
obtained from Prakash et al. jlSj. The calculation is car- 
ried out to 100 fm/c, although most of the rescattering 
finishes by about 30 fm/c. The rescattering calculation 
is described in more detail elsewhere 
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Calculations are carried out assuming initial parameter 
values and particle multiplicities for each type of parti- 
cle. In the last stage of the calculation, the freeze-out 
and decay momenta and space-times are used to pro- 
duce observables such as pion, kaon, and nucleon mul- 
tiplicities and transverse momentum and rapidity dis- 
tributions. The values of the initial parameters of the 
calculation and multiplicities are constrained to give ob- 
servables which agree with available measured hadronic 
observables. As a cross-check on this, the total kinetic 
energy from the calculation is determined and compared 
with the RHIC center of mass energy of -^/s = 130 GcV 
to see that they are in reasonable agreement. Particle 
multiplicities were estimated from the charged hadron 
multiplicity measurements of the RHIC PHOBOS exper- 
iment p^ . Calculations were carried out using isospin- 
summed events containing at freezeout about 5000 pions. 
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FIG. 1: Transverse mass distributions from the rescattering 
model. The lines are exponential fits to the distributions and 
the slope parameters are shown. 



500 kaons, and 650 nucleons. The hadronization model 
parameters used were T — 300 MeV, ay— 2 A, and Thad—^ 
fm/c. It is interesing to note that the same value of Thad 
was required in a previous rescattering calculation to suc- 
cessfully describe results from SPS Pb+Pb collisions [ pT| . 
Figure 1 shows mx distributions for pions, kaons, and nu- 
cleons from the rescattering calculation for 6 = fm near 
midrapidity (— 1 < y < 1) fitted to exponentials of the 
form exp{—mT/B), where B is the slope parameter. The 
extracted slope parameters shown in Figure 1 are close in 
value to preliminary measurements from the STAR ex- 
periment for the TT", , and anti-proton of 190 ± 10, 
300 ± 30, and 565 ± 50 MeV, respectively 0. Thus, 
we see that if all hadrons begin at a common tempera- 
ture of 300 MeV, the hadronic rescattering alone is able 
to generate enough radial flow to account for the differ- 
ences in slope among the pion, kaon, and nucleon mx 
distributions. 

The elliptic flow and two-pion HBT observables are 
also calculated from the freeze-out momenta and space- 
time positions of the particles at the end of the rescat- 
tering stage. The elliptic flow variable, V2, is defined as 




V2 = {cos{2(f))); (f) = aTcta,ii{py / px) 



(4) 



where Px and Py are the x and y components of the par- 
ticle momentum, and x is in the impact parameter direc- 
tion and y is in the "out of plane" direction (i.e. a; — z is 
the reaction plane and z is the beam direction). The HBT 
pion source parameters are extracted from the rescatter- 
ing calculation using the same method as was applied for 
previous SPS-energy rescattering calculations |l^. The 
Pratt-Bertsch "out-side-long" radius parameterization is 
used ||l^, ^ yielding the four parameters Rrside, Rxout, 
Rlohq, and A, which represent two mutually perpendic- 
ular transverse (to the beam direction) radius parame- 
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FIG. 2: Comparison of «2 calculated from the rescattering 
model for & = 8 fm with STAR measurements for pions and 
nucleons. The plotted points with error bars are the rescat- 
tering calculations and the lines show the trends of the STAR 
measurements. Average errors on the STAR measurements 
are < 0.002 for pions and 0.006 for protons-|-antiprotons. 



ters, a radius parameter along the beam direction, and 
a parameter related to the "strength" of the two-pion 
correlations, respectively. 

Figure 2 shows the pr dependence of V2 for pions and 
nucleons extracted from the 6 = 8 fm rescattering calcu- 
lation compared with the trends of the STAR measure- 
ments for TT"*" -|- 7r~ and p + p — bar at 11 — 45% centrality 
101, which roughly corresponds to this impact parameter. 
As seen, the rescattering calculation values are in rea- 
sonable agreement with the STAR measurements. Thus, 
the same rescattering mechanism that can account for 
the radial flow seen in Figure 1 also can account for the 
magnitude and pt dependence of the elliptic flow for pi- 
ons and nucleons. 

The pion source parameters extracted from HBT anal- 
yses of rescattering calculations for three different impact 
parameters, 6 = 0, 5, and 8 fm, are compared with STAR 
7T~ measurements at three centrality bins j^] in Figure 
3. Note that the PHENIX HBT results |] are in basic 
agreement with the STAR results. The STAR centrality 
bins labeled "3", "2", and "1" in the figure correspond 
to 12% of central, the next 20%, and the next 40%, re- 
spectively. These bins are roughly approximated by the 
impact parameters used in the rescattering calculations, 
i.e. the average impact parameters of the STAR central- 
ity bins are estimated to be within ±2 fm of the rescatter- 
ing calculation impact parameters used to compare with 
them. In the left panel, the centrality dependence of the 
HBT parameters is plotted for a pt bin of 0.125 — 0.225 
GeV/c. In the right panel, the tot dependence of the 
HBT parameters is plotted for centrality bin 3, for the 
STAR measurments, or 6 = fm, for the rescattering 
calculations. Although there are differences in some of 
the details, the trends of the STAR HBT measurements 



are seen to be described rather well by the rescattering 
calculation. 

As shown above, the radial and elliptic flow as well as 
the features of the HBT measurements at RHIC can be 
adequately described by the rescattering model with the 
hadronization model parameters given earlier. The re- 
sults of the calculations are found to be sensitive to the 
value of Thad used, as was studied in detail for SPS rescat- 
tering calculations Q. For calculations with Thad > 1 
fm/c the initial hadron density is smaller, fewer collisions 
occur, and the rescattering-generated flow is reduced, re- 
ducing in magnitude the radial and elliptic flow and most 
of the HBT observables. Only the HBT parameter RLong 
increases for larger Thad reflecting the increased longitu- 
dinal size of the initial hadron source, as seen in Equa- 
tion 3. One can compensate for this reduced flow in the 
other observables by introducing an ad hoc initial "flow 
velocity parameter", but the increased RLong cannot be 
compensated by this new parameter. In this sense, the 
initial hadron model used in the present calculations with 
Thad ~ 1 fm/c and no initial flow is uniquely determined 
with the help of RLong- 

At this point, it is appropriate to discuss how physical 
the initial conditions of the present rescattering calcula- 
tion are. In order to use this picture, one must assume: 
1) hadronization occurs very rapidly after the nuclei have 
passed through each other, i.e. Thad — 1 fm/c, 2) the 
hadrons thermalize rapidly, and 3) hadrons or at least 
hadron-like objects can exist in the early stage of the 
collision where there are maximum values of T and p of 
300 MeV and 8 GeV/fm^, respectively. 

Addressing assumption 3) first, in the calculation the 
maximum number density of hadrons at mid-rapidity at 
t — fm/c is 6.8 fm"^, rapidly dropping to about 1 
fm^^ at i = 4 fm/c. Since most of these hadrons are 
pions, it is useful as a comparison to estimate the effec- 
tive volume of a pion in the context of the tt — tt scat- 
tering cross section, which is about 0.8 fm? for s-waves 
p^ . The "radius" of a pion is found to be 0.25 fm and 
the effective pion volume is 0.065 fw? , the reciprocal of 
which is about 15 fm~^. From this it is seen that at 
the maximum hadron number density in the calculation, 
the particle occupancy of space is estimated to be less 
than 50%, falling rapidly with time. One could specu- 
late that this may be enough spacial separation to allow 
individual hadrons or hadron-like objects to keep their 
identities and not melt into quark matter, resulting in a 
"super-heated" semi-classical gas of hadrons at very early 
times, as assumed in the present calculation. 

Assumptions 1) and 2) can both be motivated by the 
Color Glass Condensate model Before the colli- 

sion of two relativistic nuclei, the nuclei can be described 
as two thin (lorentz-contracted) sheets of Color Glass (a 
dense glass of gluons). The use of the term "glass" is 
not merely an analogy but is mathematically rigorous 
due to the time-dilated nature of the gluon field which 
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FIG. 3: Comparison of HBT source parameters from rescattering with STAR measurements as a function of centrality bin 
(see text) and rriT- The STAR measurements are the solid symbols and the rescattering calculations are the open symbols. 
The errors on the STAR measurements are statistical+systematic. 



behaves like a liquid on long time scales but a solid on 
short ones In the usual picture, after the collision 

takes place the Color Olass melts into quarks and p-Iuotis 



in al timescale of about H fm /c at RHTC energy^ and 



then the matter expands and thermalizes into quark mat- 
ter by about 1 fm/c. The timescale in the Color Glass 
Condensate model for thermalized matter matches the 
timescale needed in the rescattering model for an initial 
thermalized hadron gas. Thus one is tempted to modify 
the collision senario such that instead of the Color Glass 
melting into quarks and gluons just after the collision, the 
sudden impact of the collision "shatters" it directly into 
hadronic fragments which then thermalize on the same 
timescale as in the parton senario due to the hadronic 
strong interactions. 

In summary, a thermal model + hadronic rescattering 
picture is able to adequately describe the large elliptic 
flow and small HBT radii recently measured at RHIC. A 
feature of this picture is a very early hadronization time 
of about 1 fm/c after the initial collision of the nuclei. 
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